Resonating color state and emergent chromodynamics in the kagome antiferromagnet 
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We argue that the spin-wave breakdown in the Heisenberg kagome antiferromagnet signals an 
instabihty of the ground state and leads, through an emergent local constraint, to a quantum 
dynamics described by a gauge theory similar to that of chromodynamics. For integer spins, we 
show that the quantum fluctuations of the gauge modes select the \/3 x \/3 Neel state with an on-site 
moment renormalized by color resonances. We find non-magnetic low-energy excitations that may 
be responsible for a deconfinement "transition" at experimentally accessible temperatures which we 
estimate. 



Spin liquids are intriguing states of matter that do 
not break any symmetry at zero temperature, and were 
conjectured by Anderson in the form of a resonating 
valence bond (RVB) liquidii Contrary to conventional 
magnets which have long-range order and spin wave 
excitations described by an effective non-linear sigma 
model,— RVB states have emergent gauge excitations.— 
Gauge degrees of freedom appear because of local con- 
straints within a low-energy manifold, such as the hard- 
core constraint for singlets^ or, more generally, some "ice- 
rules".— As a result, the magnetism can be described in 
terms of "electrodynamics" and Coulomb phases which 
find striking experimental realizations in spin-ice systems 
at finite temperatures.— Recently, experiments on the 
two-dimensional spin- 1/2 Heisenberg kagome compound 
ZnCu3(OH)6Cl3 suggest a possible quantum spin-liquid 
state at zero tcmpcrature<^ Theoretically for S = 1/2 
systems, there is indeed evidence for a spin-liquid state,- 
not inconsistcnt^i^ with a gapless RVB state.~ From a 
different route, the large-S approach provides important 
insights, such as integer vs. half-integer spin effectsJ^ 
For the kagome system, the large-S" spin-wave theory 
is known to break down because of local (weathervane) 
modes with diverging fluctuationsiii"— This is the con- 
sequence of a macroscopic number of classically degen- 
erate states that can be seen as the 3-colorings of the 
kagome lattice with spins pointing at 120° apart (Fig. 
1). Contrary to spin- ice models where the constraint 
is enforced by strong local anisotropics;^ here the local 
(color) constraint arises dynamically because of an order- 
by-disorder mechanism.—"— In this paper, we show that 
the spin-wave breakdown is the signature of emergent 
gauge degrees of freedom that arc governed by a Hamil- 
tonian similar to that of quantum chromodynamics. The 
issue is whether the present theory stabilizes a resonating 
color "crystal" (with a finite Neel order-parameter and 
color resonances) or leads to a delocalization in the low- 
energy manifold, thus realizing a resonating color "spin- 
liquid" (instead of RVB). 

First, we derive the effective lattice gauge theory in the 
low-energy manifold of 3-colorings starting from largc-S*. 
It is known both from classical Monte Carlo simulations 
{S = oo)i^'i^ and from exact spectra at finite that 
the 3-coloring states (spin-ice subspace) are relevant at 
low energy. To remain in the low-energy manifold, single 



FIG. 1: (color online.) Loop-induced low-energy dynam- 
ics within the 3-coloring manifold of the kagome lattice 
(spin/color equivalence is shown), described as "gluon" dy- 
namics. Right: the q\q\ ("meson") and qiq2q3 ("baryon") 
are high-energy "matter" excitations violating the local con- 
straint. 



spin-flips are not allowed and the spin motion consists 
of the tunneling of a collective loop of spins across an 
energy barrier separating degenerate statesJ^ The loops 
are closed loops of L spins of two colors (Fig. 1). Classi- 
cally, they can be rotated freely about the effective field 
of the third-color spins, for any arbitrary collective angle 
(j). However, this local rotation can be viewed as a de- 
fect for the spin-waves and, therefore, increases the zero- 
point energy. In the Born-Oppcnheimcr approximation, 
this generates an effective energy barrier for the slow loop 
motionJ^iii We now calculate the energy barrier starting 
with the quantum Heisenberg Hamiltonian 



H — J ^ Si ■ Sj, 



(1) 



<h3> 



with J the coupling between nearest-neighbor spins S. 
To do so, we start with a classical ground state and con- 
sider a given loop (see Fig. 1 for example) rotated by 
(j). The Holstein-Primakov spin-wave equations of mo- 
tion are set up in real space (given that translation in- 
variance is broken) and diagonalized for different (j) (up 
to = 972 sites). Wc compute the zero-point energy 
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where I runs over all spin- wave energies, ui . It is a double- 
well potential with two minima located at ±7r/2 (Fig. [5]) 
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FIG. 2: (color online.) Left: Double-well potential for the 
loop motion and quantized levels with even (solid line) and 
odd (dashed line) wavefunctions. Right: Barrier height distri- 
bution as a function of loop length in the 3-coloring manifold. 
The curves shown are simple cosine functions. 

corresponding indeed to the selection of 3-colorings. The 
potential is well fitted by an inverted parabolic function, 
which we will use in the following. Repeating the cal- 
culation for loops of arbitrary lengths in a Monte Carlo 
sampling of the manifold, we have systematically com- 
puted the barrier heights Vq (Fig. [2]). For small loops, 
we find Vq = rjJSL with rj = 0.14, in agreement with 
Ref. [H; the barriers are rather small. For larger loops, a 
dispersion about the mean value is observed, which tells 
us that the barrier heights depend on the configuration. 
Note that the symmetry about L = N/3 occurs because 
rotating a loop of L = 2N/3 results in choosing another 
global plane at no energy cost. We also point out that 
anharmonic corrections will give additional contributions 
~ J5'2/3 _i6^i8,i9 rp]^^ kinetic energy of a loop is obtained 
by considering the spins of the loop in the effective field 
of the fixed third color spins (say along ~z). The Gold- 
stone mode corresponds to the collective rotation about 
z and is taken scparatly»^ Together with Vl we obtair^^i 
the hamiltonian for a single loop 

HL = ^{S'-Ay + VL{cj^) (3) 

where = X^iLi = —id/d(j> is the total spin of the 
loop, conjugate to 0, and x = L/{4:J) the moment of 
inertia. Hl is similar to an Aharonov-Bohm ring in a 
double- well potential. The flux A = LS/2 is the classical 
magnetization of the loop. For the kagome lattice, the 
only possible loops have L = 4n + 2 (n an integer) so 
that A = {2n+l)S. By considering the Berry phase 27r^, 
von Delft and Henley predicted a destructive interference 
for half- integer SJ^ Indeed, the first term of ^ is min- 
imized for = A if A is an integer or = A ± 1/2 
if A is half- integer ("Kramers" doublet). In the pres- 
ence of the external potential, this remains true. By us- 
ing a simple gauge transformation 'i'{4') = e~'''^'^u{(j)), 



u((j)) satisfies the Schrodinger equation with either peri- 
odic (integer S) or antiperiodic (half-integer S) bound- 
ary conditions. Vl{4>) has a periodicity of tt so we ex- 
pect Bloch states with energies Ek- Periodic (resp. anti- 
periodic) boundary conditions selects the bonding k = 
and anti-bonding fc = 1 states (resp. k ~ ±1/2), which 
have energy difference t^. Given that Vl{4>) = VL{~(t>), 
we have E-k = Ek and the two states k = ±1/2 remains 
degenerate [t^ ~ 0). 

To extract t^^ we solve exactly the Schrodinger equa- 
tion in the double- well parabolic potential with parabolic 
cylinder wavefunctions (Fig. [5]). In the semi-classical 
limit (large S or large L), the solution coincides with the 
Wcntzel-Kramers-Brillouin approximation and we obtain 

iL = aJexp(-^V^^L) (integers) 
(half-integer S) 

where a = Q.2Ql{i^Sf^'^^L^/^ and taking = Vo/JSL, 
i.e. neglecting the dispersion of the barrier heigths. 
Quantum tunneling of loops of size L > Lc = 
4/7r^-\/?7S'/2 is therefore suprcssed. 

We now consider the resulting low-energy dynamics on 
the lattice. In this derivation, the lowest energy states 
are the 3-coloring quantum coherent states |C).— The 
number of states scales like W ^ 1.13^ in the thermo- 
dynamic limit^i^ and is much smaller than that of the 
original spin problem ([1]), (25* +1)''^. The dynamics of 
the loops is described by the quantum Hamiltonian 

He = -J2 tc,c'\C){C'\+J2Ec\C){C\ (5) 

C,C' c 

where tc.c" ~ if the two states \C) and |C") arc con- 
nected by the tunneling of a loop of size L and tc.c = 
otherwise. The model depends only on 5" and the largest 
energy scale is = t. Anharmonicit y ^ ^i ^ ^ and further 
neighbor couplings generate a finite Ec and phase com- 
petitions that are beyond the scope of the paper— (here 
Ec = 0) . To rewrite Hf. exactly in terms of exchange op- 
erators of three colors, we have to use the 3x3 Gell-Mann 
matrices A^;^ 2 3 °f SU(3) algebra. We find 
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He = E E K^K,2K3KaK,K,, + h.C. + • • • (6) 

O c=l 

where the dots stand for longer loops. The Hamilto- 
nian (|6]) is a pure lattice gauge model describing the 
loop dynamics, and is similar to that of quantum chro- 
modynamics with the local SU(3) symmetry broken to 
U(l)xU(l).— The weathervane modes can therefore be 
seen as "gluons" . Note that the use of the eight operators 
of SU(3) is necessary because here, contrary to Ref. [2^, 
the Hamiltonian respects the global Z3 symmetry and 
flips the three types of loops. If we pursue the anal- 
ogy, the finite-energy excitations are defects in the ice 
rule that consist of (i) flipping two colors along a finite 
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stringj22, it can be seen as a "quark-anti-quark" qiqi (me- 
son) (Fig[T]). It terms of spin, it is a non- magnetic pair of 
domain- walls with irrational magnetization •\/3(+ei, — e^) 
[e-i are three unit- vectors with ei -I- 62 -I- 63 = 0). (ii) per- 
muting a triangle, this generates three quarks in a color 
singlet qiq2(l3 or 919293 (baryon), or, in the spin language, 
a non-magnetic excitation ■\/3(ei, 62, 63). These are ex- 
actly the first excitations of the Potts model, whereas 
there are other excitations for the Heisenberg model. We 
will address below the issue of confinement of these ex- 
citations. Finally, we point out that has a complex 
topological structure separated in Kempe^^ and topolog- 
ical sectorsi^ While topological sectors can be connected 
only by hipping winding loops, Kempe sectors are ensem- 
bles of topological sectors that cannot be connected by 
such global moves i^i^S 
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FIG. 3: (color online.) Finite-size scaling of the order- 
parameter of the \/3 X \/3 Neel state, extrapolating to a finite 
value (exact diagonalization and quantum Monte Carlo). In- 
set: total energy and fits assuming gapped (solid lines) or 
gapless (dashed lines) excitations. 

The ground state is in general a resonating color state, 
^c\C) in a given sector and the issue is whether 
it breaks the Z3 symmetry or not. There is an exten- 
sion of the model for which is known exactly. By 
adding Ec = ULnL{C) which penalizes states with a 
large number nL{C) of loops of size i, one obtains a 
Rokhsar-Kivelson (RK) model for color dimers<^ In this 
case, the ground state is the resonating color "spin- 
liquid", 4* = X^c l*-^) each topological sector at the 
RK point, Ul = We study the model ^ with 

Ec = by Lanczos exact diagonalization of clusters of 
size TV = 27,36,54,81,108 and T = quantum Monte 
Carlo (up to = 675), since He is free of the sign prob- 
lem. We find that (i) "F belongs to the largest topologi- 
cal and Kempe sector (which is non degenerate and has 
all topological numbers equal to //3 where I is the lin- 
ear dimension), (ii) 4' has larger amplitudes onto the 
six X -^3 states, and a finite associated Neel order- 
parameter, rn^ which extrapolates to a finite value for 
infinite size, 



(Fig. [3]). To confirm the Neel order, we also calculate 
the low-energy spectrum by exact diagonalization. For 
such a discrete broken symmetry, we expect six quasi- 
degenerate states separated from higher energy states. 
While there is no clear separation of scale, the gap does 
decrease as exp(— aA^^^^) which reflects the development 
of a macroscopic energy barrier (Fig. |4|) . The low-energy 
spectrum is therefore perfectly compatible with a Neel or- 
der with the same pattern as predicted earlie r^^i^^i'^^ and 
on-site magnetization renormalized by color fluctuations. 
However, there are unconventional non-magnetic excita- 
tions at low-energy of order t (gluons) which bear some 
similarity with the low-lying states of exact spectraJ^i^^ 
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FIG. 4: (color onUne.) Left: energy spectrum (A'^ = 108) 
showing no clear separation of scale. Right: the lowest finite- 
size gap follows AE/t ~ exp(— aA*''^'^^) confirming the broken 
symmetry of the \/3 x \/3 Neel order. 

To discuss the consequences of these unconventional 
excitations and analyze these results more quantitatively, 
we consider the effective gauge theory in the continuum 
limit. The ice-rule provides an exact mapping of the 
discrete states onto two-component height vectors h(r) 
living on the dual latticej^i^i^ For instance the six 
V3 X states map onto flat interfaces. At T = 0, 
the simplest action compatible with the symmetries is a 
sine-Gordon model, 



S = 



K 
2? 



(a,h)2 + c2(vh)= 



3 

a=l 



cos(Qa.h) 

(7) 

where K is the stiffness, the amplitude of the lock- 
ing potential, and |Qq-| ~ An/y/SM^^^^ Classical min- 
ima of ^ correspond to the "flat" -v/S x a/3 Neel states, 
with zero-point oscillations described by the two "gluon" 
modes Wq = (c^q^ + A2)V2. At T = and in (2+1) 
dimensions, no quantum phase transition is expectedi^ 
The gluons mediate an interaction between the defects 
(quarks) that is conflning if A is flnite as V{r) ~ Ari^ 
To obtain these parameters, we have used the flnite-size 



0.63 — 0.68; depending weakly on Lc correction to the total energy. 



N 



eo - 0.6586 



2c 
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a hexagonal lattice of A'' sites, without a gap (dashed 
line in Fig. 2]) . We clearly see that the fit is improved at 
large distance (for a finite Lc) by using j^l^2 exp(— //^) 
where f = c/A gives the confinement scale (solid line), 
and c l3At, ^ = 15.5, eo = -0.1430t. To es- 
timate the stiffness we calculate the order-parameter, 
moo = (e''^ '') with |Q| 4^/3,^2,^^ Assuming a small 

gap, moo = e^^*^ ^'^ ^ = and by using rrioo ^ 0.65, 

we extract K ^ 9Qt. 

Thermal effects. We first discuss the thermal restora- 
tion of the discrete symmetry of the effective model: the 
interface described by ([7]) undergoes a roughening tran- 
sition in the Kosterlitz-Thouless universality class (trig- 
gered by the condensation of gluons) at a critical tem- 
perature usually obtained from the scaling dimension of 
the locking potential, 2 — IQal^T/AirK , leading to Tc = 
3K/{2tt). For T < Tc the system breaks the discrete Z3 
symmetry with a V3 x -s/S state, the weathervane modes 
(gluons) are gapped and the defects (quarks) are con- 
fined. For T > Tc the Z3 symmetry is restored, it is a crit- 
ical quark-gluon plasma with gapless gluons (77 = 4/3j^ 
hence a dynamical response x"{'^) ~ o;''/^"^). For in- 
teger spins, we estimate Tc/J — 0.16 (S ~ 1), 0.06 
(S' = 2), 0.02 {S 3), etc. by using K - 90t and 
Eq. m For half-integer spins, since t = (at the one-loop 
order) and t t^ (with two- loop co-tunneling)^ the 



system is immediately in its infinite temperature phase 
and is critical. Could this deconfinement transition man- 
ifest itself experimentally? The transition should indeed 
occur in the corresponding spin-ice system with a dis- 
crete symmetry. In the model ([1]) with SU(2) symmetry, 
the Mermin- Wagner theorem ensures that no continuous 
symmetry breaking will occur at finite temperature: the 
propagating spin- waves induce a finite correlation length, 
£.2(1 ^ cxp{p/T)^ which cuts off the infrared divergence 
of the Kosterlitz-Thouless transitioni^ Strictly speaking 
Tc is a crossover temperature between two paramagnets. 
However, the usual thermodynamic singularities will re- 
main sharp because ^ exp(JS'^/rc) is very large at 
T = Tc. 

In conclusion, we argued that the Heisenberg kagome 
antiferromagnet has long-range Neel order at T = 
for integer spins. Its magnetism is unconventional due 
to an underlying emergent chromodynamics: there are 
color fluctuations and additional low-energy excitations 
described as gapped "gluons" which may condense at Tc 
("quark" deconfinement). For half-integer spins, Tc = 
at the lowest order. The difference between half-integer 
and integer spins refiects the essential quantum nature of 
the local modes. 
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